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High-energy isolated attosecond pulses required for the most intriguing nonlinear attosecond 
experiments as well as for attosecond-pump/attosecond-probe spectroscopy are still lacking 
at present. Here we propose and demonstrate a robust generation method of intense isolated 
attosecond pulses, which enable us to perform a nonlinear attosecond optics experiment. By 
combining a two-colour field synthesis and an energy-scaling method of high-order harmonic 
generation, the maximum pulse energy of the isolated attosecond pulse reaches as high as 
1.3 jiJ. The generated pulse with a duration of 500 as, as characterized by a nonlinear auto- 
correlation measurement, is the shortest and highest-energy pulse ever with the ability to 
induce nonlinear phenomena. The peak power of our tabletop light source reaches 2.6 GW, 
which even surpasses that of an extreme-ultraviolet free-electron laser. 
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Attosecond science has emerged as an important research 
area of ultrafast phenomena within the past decade^ 
Owing to its numerous successes, attosecond science has 
created important new knowledge (see review ref. 2) of the 
fundamental science of the interaction between electrons and 
photons. Recently, some key research topics have been identified, 
which are expected to make major breakthroughs for the next 
attosecond frontiers, such as the MHz high-repetition high-order 
harmonic (HH) frequency combs for the joint frontier of 
precision spectroscopy and ultrafast science, attosecond-pump/ 
attosecond-probe experiments for observing and controlling 
electronic processes in atomic and molecular physics, and 
nonlinear science occurring on the attosecond time scale. 

To seriously tackle the above interesting research topics, one of 
the most important issues is the development of high-power 
isolated attosecond pulses (lAPs) and/or attosecond pulse trains. 
lAPs and attosecond pulse trains are produced using HH 
generation (HHG)^'"* in gases. To date, high-power attosecond 
pulse trains have successfully been generated owing to research 
on HH energy scaling using a loose-focusing geometry^. 
Attosecond pulse trains with sufficiently high energy in the 
microjoule range with a high conversion efficiency on the 10""* 
level were generated^, which are intense enough for implementing 
some applications, such as nonlinear attosecond optics^'^ without 
the assistance of laser pulses, single-shot femtosecond 
holography^ and an external injector for a seeded free- 
electron laser (PEL) in the extreme-ultraviolet (XUV) region. In 
contrast, the output energy of lAPs is still not sufficient for 
various applications, although the shortest pulse duration 
achieved is sub- 100 as in ref. 11. The widespread application of 
lAPs has been limited because of the low photon flux and the 
complexity of the laser systems required to produce lAPs (see 
review in ref. 12). At present, the highest lAP energy is less than a 
few tens of nanojoules^^ at a 27-eV photon energy. The highest 
conversion efficiency demonstrated is only 10~^to 10~^, even if 
Xe gas is used. Consequently, the production, characterization 
and application of high- energy lAPs are still under active 
investigation. 

Here we report a harmonic generation method suitable for 
nonlinear pump/probe experiments that demonstrates the 
capability of time-resolving attosecond dynamical processes with 
lAPs. Our generation scheme is robust and straightforward for 
scaling up the HH energy, which is based on an infrared two- 
colour (TC) laser field synthesis^^ and the energy- scaling 
method^. By carefully designing the HHG configuration, we can 
obtain lAPs with a record energy of up to 1.3 |ij at ~ 30 eV, thus 
demonstrating a > 100-fold energy enhancement compared with 
previous reports The conversion efficiency achieved is 
improved to 1.1 x 10""*, owing to favourable phase-matching 
technique. Even though our TC method leads to weak multiple 
burst emissions occurring for specific combinations of the 
carrier-envelope phase (CEP) and relative phase, we nevertheless 
can observe the nonlinear signal corresponding to lAPs, because 
the nonlinear process automatically selects predominantly the 
contributions from the most intense lAPs generated for certain 
favourable phase combinations. In fact, we present the direct 
temporal characterization of an lAP by a second-order 
autocorrelation (AC) measurement using nonlinear phenomena 
in N2 molecules. The AC measurement is one of the simplest 
applications of a nonlinear pump/probe experiment. The 
measured AC traces indicate not only the temporal duration of 
the interacting lAPs but also provide evidence for the nonlinear 
interaction with the target medium. From the 500-as pulse 
duration determined by the AC measurement, the lAP's peak 
power can be evaluated to be 2.6 GW with weak satellite pulses at 
the generation point. Furthermore, in our TC-HHG experiment. 



we create a mid-plateau region exhibiting a quasi-continuous 
spectrum with an energy of 10 jij. From the AC measurement of 
this quasi- continuous spectrum at the mid-plateau, we evaluate 
that a quasi-IAP with pulse duration of 375 as with ~0.27 
intensity satellite pulses at ± 6.7 fs dominated the nonlinear 
interaction with the N2 target. 



Results 

Harmonic spectra and output energy. We strictly considered 
and designed the configuration of a thoroughly engineered gen- 
eration scheme (see Methods section) for scaling up the HH 
output energy. Figure 1 shows typical single- shot HH spectra of 
Xe gas as function of photon energy. The focused intensity in the 
TC case (/tc) was ~5 x lO^^Wcm"^ with an intensity ratio of 
(^~0.15. The HH yield is optimized by adjusting the focusing 
point with respect to the gas cell and by varying the gas pressure. 
As the difference between the dispersions for the 800- and 1,300- 
nm pulses in Xe is negligible, the phase-matching condition can 
be satisfied by optimizing the gas pressure. The HH yield gra- 
dually increases as the Xe gas pressure increases, and exhibits a 
quadratic dependence on pressure. This means the conversion 
efficiency and the photon flux of HHG depend on the Xe gas 
pressure. The cutoff intensity of TC-HHG is optimized at a gas 
pressure of ~0.8Torr, which is in good agreement with the 
optimum gas pressure predicted by the phase-matching theory^ ^ 
under these experimental conditions. 

When the one-colour (OC) field is focused to an intensity of 
7 X 10^^ W cm ~ ^, the HH spectrum has a discrete structure (blue 
curve). In the TC scheme (red curve), a quasi-continuous HH 
spectrum appears in the lower- order region, whereas the quasi- 
continuous HH disappears at harmonic orders higher than the 
17th. We can clearly see a smooth continuum HH spectrum 
around the HH cutoff region. Note that the intensity of the 
TC-HH is comparable to that of the OC-HH. From the 
wavelength-scaling law of the TC-HHG (see Fig. 8) under the 
= 0.15 condition, the conversion efficiency decreases only by a 
factor of (1,300/800) -^-^-0.71 compared with the OC-HHG 
case. The maximum HH photon energy is slightly extended to 
35 eV even though the focused intensity is lower than for the 
OC-HHG. This photon energy extension is explained by the 
cutoff formula for TC-HHG^^. 

HH order for 800 nm 
9 11 13 15 17 19 21 23 25 
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Photon energy (eV) 

Figure 1 | Experimentally obtained single-shot HH spectra from Xe. Blue 
curve: OC laser field (800 nm: 11 m J); red curve: TC laser field (800 nnn: 
9 nnJ + 1,300 nm: 2.5 mJ). Green dotted curve: reflectivity of the Sc/Si 
multilayer mirror; black dotted curve: reflectivity of the SiC bulk mirror. 
Inset: (a) measured intensity variation of the HH cutoff (~29eV: 19th HH 
order) over 200 laser shots, (b) The spatial profile of the cutoff harmonics 
of the TC-HHG. 
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The inset (a) of Fig. 1 shows the measured intensity variation of 
the HH cutoff ( - 29 eV: 19th HH order) over 200 laser shots. The 
intensity of the OC-HH has a high stability, as shown by the blue 
points, whereas the intensity of the TC-HH (red points) fluctuates 
between 1 and 0.3. As already reported in ref. 14, these variations 
are due to the intensity changes of the TC field caused by random 
fluctuations in the CEP value of the non-CEP -stabilized 
Ti:sapphire laser and slight jitter of the relative phase 
introduced by the Mach-Zehnder interferometer used for 
creating the TC fields. According to the discussion of the 
variations of these two kinds of phase effects (see Discussion 
section), the HH yield at the cutoff should be maximized, when 
the continuous spectrum frequently appears for certain 'good' 
combinations of the CEP and relative phase, whereas only weak 
HH exhibiting the quasi- continuum spectrum is generated for 
'bad' phase combinations. We have actually observed this 
correlation between the HH yield and the spectral shape in the 
experiment. Note that harmonic shots with the quasi- continuum 
spectrum hardly contribute to form the AC trace, because their 
intensity is quite low. 

The output beam divergence of the TC-HH at the cutoff was 
measured to be ^O.Smrad full-width half maximum with a 
Gaussian-like profile (see Fig. lb inset). The almost perfect 
Gaussian profile of the HH suggests that there is no density 
disturbance due to high ionization in the harmonic medium. The 
pulse energy of the continuous HH spectrum (28-35 eV) is 
evaluated to be 1.3 |ij with a conversion efficiency of 1.1 x 10""* 
from the input TC energy (800 nm: 9mJ+l,300nm: 2.5 mj). 
This value is almost 100-fold higher than the highest pulse 
energies ever reported before The generation efficiency from 
the Ti:sapphire laser is 3.7 x 10 because 30 mJ of Ti: sapphire 
energy was used for creating the TC field. Moreover, the mid- 
plateau region (14-29.5 eV) exhibiting a quasi- continuous 
spectrum in the TC-HHG has an output energy of ~ 10 jij. 

Nonlinear AC measurement for attosecond pulses. To evaluate 
the temporal characteristics of the microjoule HH fields, we uti- 
lized an AC method using a spatially split autocorrelator. We 
employed two concave mirrors (i^ = 200mm) for focusing the 



HH: one is a bulk SiC mirror for the mid-plateau and the other is 
a Sc/Si multilayer mirror for the cutoff. The reflectivity for each 
concave mirror is shown in Fig. 1 . We have adopted the N ^ ion 
yield from fragmenting N2 molecules as the indicator of the 
second-order AC signal (see Methods section). This choice is 
motivated by the findings of our earlier experiment^, in which 
the N ^ ion yield was sufficiently high to obtain the AC trace of 
an attosecond pulse in the relevant HH photon energy region. 
The dotted curves in the top panels of Fig. 2a,b are the measured 
AC traces of the mid-plateau TC-HH and OC-HH fields, 
respectively. The signals at each delay time are obtained by 
accumulating signals for 1 x 10^ laser shots. We can see distinct 
peaks emerge at every half period (1.33 fs) of the Ti:sapphire 
laser field in the AC trace depicted in the top panel of Fig. 2b. 
This result ensures that our AC measurement can certainly 
reveal the temporal shape in the attosecond regime. In contrast, 
the four peaks corresponding to multiple bursts on both side of 
the main peak (Ofs) are well suppressed in the TC-HH, even 
though we can observe the side peaks at ~ ± 6.7 fs (indicated by 
arrows in Fig. 2a) with a relative AC amplitude of ~0.45 
compared with the main peak. These two side peaks are 
consistent with our simulated result (see Discussion section) of 
the TC-HHG. 

To observe the pulse duration more precisely, we increased the 
time resolution of the AC measurement by reducing the delay 
step from 148 to 28 as. The bottom panels in Fig. 2 show the 
measured AC signal for OC-HH (blue points) and TC-HH (red 
points). In the mid-plateau of TC-HH, the full- width half 
maximum of the AC was measured to be Tac = (530 ± 75) as by 
fitting the trace with a Gaussian function (not shown). The pulse 
duration is thus Thh = t^ac/V^ ~ 375 as. On the other hand, 
the measured AC duration for the OC-HHG case was 
'^Ac — (490 ± 60) as. Thus, the duration of each peak in the 
attosecond pulse train is evaluated to be THH^350as, which is 
slightly shorter than the pulse duration of the TC-HH. The 
broadening of the AC duration of the TC-HH compared with the 
OC-HH case is owing to the fact that the AC trace of the TC-HH 
originates from the temporal profiles averaged over the random 
CEP and slightly fluctuating relative phase. In addition, the HH 
spectral shapes are different in both cases. 



a b 




At {is) At (is) 

Figure 2 | Measured AC traces from the side peak of N+ ion signals, (a) Mid- plateau of TC-HH; (b) mid-plateau of OC-HH. The time resolutions 
of the top and bottom panels correspond to 148 and 28 as, respectively. The error bars show the s.d. of each data point. The grey solid profiles are the AC 
traces of the simulated HH spectrum. Inset of (a): measured TOF spectrum of ions at around m/z = 14. The inverted triangles and black-filled 
inverted triangle correspond to the KE of 3 and OeV, respectively. 
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As we successfully obtained the AC trace of the TC-HH in the 
mid-plateau, we furthermore explored the temporal characteriza- 
tion for the HH cutoff. In this experiment, we switched the 
focusing mirror from the SiC mirror to the Sc/Si multilayer 
mirror. The signal is scanned every 148 as from —8.5 to 8.5 fs 
with 1 X 10^ laser shot accumulation at each delay point. The 
resultant signal points versus delay are shown in the top panel in 
Fig. 3. Three distinct signal points around 0 delay seem to be a 
correlation peak and there are no noticeable side peaks in this 
figure, although our simulation result (see Discussion section) 
exhibits side peaks at ± 6.7 fs with ~0.2 relative AC amplitude 
compared with the maximum peak at 0 delay. To confirm that the 
signal around 0 delay is certainly the correlation signal and that 
there are no significant side peaks at ± 6.7 fs, we performed a 
high -resolution AC measurement in the delay scanning ranges 
around Ofs and ± 6.7 fs, as shown in the bottom panels in Fig. 3. 
We accumulated 1x10^ laser shots at each delay point and set the 
scanning delay step At to be 28 as. We can demonstrate a clear 
peak at 0 delay as an AC in the middle-bottom panel in Fig. 3. The 
AC duration was evaluated to be Tac — (^00 ±90) as by fitting the 
experimental trace to a Gaussian function; hence, the pulse 
duration is estimated to be Thh ^ 500 as. We also find that there 
are no noticeable side peaks in the measured AC traces in the left- 
bottom and right-bottom panels in Fig. 3. This is owing to the fact 
that the background noise amplitude (see Methods section) is 
comparable to the amplitude of the side peaks expected from our 
simulation. The noise magnitude against the peak AC amplitude is 
estimated to be ~0.24; thus, we cannot clearly find side peaks 
with amplitudes <0.24. In other words, we can safely conclude 
that the side-peak ratio should be <0.24. As discussed below, we 
can find both pre- and post -pulse in the calculated temporal pro- 
file at any CEP and relative phase. By assuming that the HH pulse 
consists of three (pre-, main and post-) pulses and the upper limit 
of the side peaks in the AC trace should be 0.24, we can estimate 
the possible largest peak height ratio of the pre- and post-pulses to 
be ~0.14. Further detailed analysis of the pre-/post-pulse issue is 
presented in the next section with the help of our numerical model 
of the TC-HH and statistical analysis of the relative phase. 

From the measured pulse energy and pulse duration, the peak 
power of the lAP nonlinearly interacting with the target medium 



is evaluated to be 2.6 GW with weak satellite pulses (see 
Discussion section). This ultrahigh power from our tabletop 
light source even surpasses that from a compact XUV FEL^^' . 
The peak brightness is also estimated to be ~2 x 10^^ photons 
per s~^ mm~^mrad~^, assuming a 60-jim beam diameter at 
the exit of the gas cell. 

Discussion 

Finally, we examine attosecond pulse generation in our TC-HH 
scheme by numerical simulations. In these simulations, the 
single-atom response is calculated within the strong-field 
approximation. The propagation of the TC laser pulse and the 
high harmonics are computed separately in cylindrical coordi- 
nates All the parameters, such as laser intensity, target gas 
pressure and medium length, are determined by the experimental 
conditions. Here we assume that the CEP 0ce fluctuates 
randomly during the AC measurement, whereas the statistical 
fluctuations of the relative phase 0o are characterized by the 
measured histogram of the relative phase in the TC inter- 
ferometer. To evaluate the relative phase variation with time in 
the TC interferometer, we injected a continuous-wave (CW) laser 
beam with a wavelength of 547 nm into the TC interferometer in 
this measurement. We recorded spatial interference fringes on the 
superposed beam profile of the CW laser with a charge- coupled 
device camera (exposure time: ~ 10 jis with 10 Hz repetition rate, 
which is synchronized with the pump laser) during ~35min, 
which is equivalent to the data acquisition time of one AC trace in 
our measurements. The phase of the spatial fringes at each 
recorded profile was converted to the relative delay and relative 
phase corresponding to the 1,300-nm field. The measured time 
evolution of the relative phase, 0o> is shown in Fig. 4a. We can 
then evaluate the s.d. of 00 fluctuations, a, to be 0.2371 from this 
measurement, whereas the offset of 0o is arbitrary. We adjust the 
offset in Fig. 4a such that the average of 0o is equal to 0. To 
improve the S/N ratio of the AC trace, we overlapped ten AC 
traces in Figs 2 and 3. The a- value with a typical normal 
distribution of the ten AC traces' (po fluctuations was almost the 
same during the experiment, even though the phase profile as a 
function of time was different. 
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The histogram of our measured relative phase exhibits a typical 
normal distribution expressed by ^ exp[— (0q — (0Q))^/(2cr^)] 
with a s.d., cr, of 0.2371 (see Fig. 4b). The mean value of (/)o> denoted 
as <0o> in this equation, cannot be determined by our jitter 
measurement. Nevertheless, we set <(/)o> to be 0.1471 rad in our 
averaging procedure in the calculations for reasons discussed 
below. Using the parameters mentioned above, we first calculate 
the temporal intensity profile for each 0o (calculation step: 
lOOmrad) value, varying 0ce in the range from — 7i/2 to 7i/2; 
later, we average the calculated HH intensity profile. We call the 
temporal profile of this averaged field the 'CEP -averaged temporal 
profile' at each 0o (see Fig. 5). Note that we adjust the offset of the 
time axis of each temporal profile of the HH field so as to locate 
the highest peak at the time origin before CEP averaging. This 
adjustment is consistent with the temporal profile to be measured 
with the spatially split AC technique, because the delay time {At) 
at the highest AC signal does not indicate the absolute time of the 
HH emission with respect to the driving TC-HH field, but 
coincides with the zero delay time between two replicas of the 
measured pulse. In an actual setup (see Fig. 9), we produce two 
time- delayed replica pulses from one HH pulse using two silicon 
beam splitters placed at Brewster's angle for the driving field^^. 



After performing the CEP average, we superpose the 'CEP- 
averaged temporal profiles' at each 00 (see Fig. 5) with the above- 
mentioned weighting factor of the normal distribution of 0o (see 
solid black curve in Fig. 4b). We have numerically searched for 
the mean value, <0o>, in the normal distribution so as to 
minimize the side peaks of the superposed temporal profile, and 
then found <(/)o> to be 0.1471 rad as mentioned above. This 
numerical procedure is compatible with our experimental app- 
roach of adjusting the relative phase in the TC interferometer 
such that the TC-HH exhibits the continuous spectrum most 
frequently. Thus, it is justified to set <0o) to be 0.1471 rad in the 
normal distribution. The resulting 0o-weighted superposition of 
the 'CEP-averaged temporal profile', which we call the 'averaged 
temporal profile', in the mid-plateau region is shown as red-filled 
curve in Fig. 6a and that in the cutoff region is shown as green- 
filled curve in Fig. 6b. The averaged temporal profile includes 
both the effect of random fluctuation of 0ce and the measured 
relative phase jitter of 0o- The overall selected bandwidths of the 
mid-plateau and cutoff harmonics are 13-29 eV and 27-36 eV, 
which are determined by the reflectivity of the SiC and Sc/Si 
mirrors, respectively. Note that 13-23 eV HH effectively con- 
tribute to creating the temporal profile in the mid-plateau 
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Figure 4 | Evaluation of the relative phase, (a) Measured relative pliase versus time, (b) A histogram obtained from a with a normal distribution fit 
(black solid curve). The intensity ratio of the CEP-averaged profile of the main (blue solid curve), pre- (red dotted curve) and post-pulse (red solid curve). 
The grey-filled area corresponds to the ±cr area. 
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TC-HH profile due to the reflectivity edge of the SiC mirror. We 
notice that two satellite pulses appear at ± 6.7 fs in both TC-HH 
profiles. The relative intensity of the higher satellite pulse 
compared with the main pulse for the mid-plateau TC-HH is 
0.27. In contrast, the relative intensity of the cutoff TC-HH is 
considerably decreased to 0.1. Note that this side-peak intensity 
can be decreased by extracting slightly higher cutoff 
harmonic orders (>29eV). We also observe a background 
pedestal under the main pulse in both TC-HH profiles. This small 
pedestal originates from weak multiburst emission around the 
centre of the TC electric field at the 'wrong' combinations of CEP 
and relative phase. The intensity of this pedestal is very low, 
because the HH yield is also low. Even though the temporal 
profile has side peaks ( ~ 0.27 relative intensity) at ± 6.7 fs, the 
mid-plateau of TC-HH results in an lAP profile having a 390-as 
temporal duration for the central pulse. When we extract the 
cutoff harmonics, the expected temporal duration is 525 as 
with weak satellite pulses. We also calculated the AC trace by 
averaging over the CEP and relative phase according to the 
procedure shown in Fig. 5, the simulated results are shown by 
dashed curves in Fig. 6 and solid curves in Figs 2 and 3. As shown 
in these figures, the calculated AC traces are in good agreement 
with the experimental data. We can clearly find a significant 
reduction of the pulses forming an attosecond pulse train with a 
period of 1.33 fs in both TC-HH fields by comparing these 
profiles with that of the OC-HH in the mid-plateau region. Thus, 
we can conclude that the intensity of multiple burst emissions in 
our TC-HHG can be significantly reduced compared with another 
gating method^^'^^ In fact, the previously reported AC profile^ ^ 
exhibits a duration in the femtosecond region when a non-CEP- 
stabilized multicycle driving laser was utilized for generating 
lAPs. We notice that the intensity ratio between the main (0 fs) 
and the satellite pulses ( ± 6.7 fs) of the AC traces is much higher 
(approximately twice) than that of the temporal profiles shown in 
Fig. 6, respectively. As the temporal profile of the TC-HH cutoff is 
composed by three pulses, the simulated AC intensity ratio of the 
side peak is ~ 0.2 compared with the main pulse, even though the 
temporal intensity ratio of the side peak is ~0.1. 

In addition, we discuss in detail how these pre- and post- 
pulses at ± 6.7 fs affect the measured temporal profiles with the 
statistical fluctuations of these phases. Here we do not consider 
the background pedestal under the central main pulse (in the 
region within ±2.7fs) originating from weak multiple burst 
emissions. In Fig. 4b, the peak intensity of the main pulse, the 
pre- pulse and the post-pulse are plotted versus 00 as blue solid, 
red dotted and red solid curves, respectively. We can see that the 
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Figure 6 | Averaged temporal profiles and their AC traces in TC-HHG. 

(a) Mid-plateau of TC-HH; (b) cutoff of TC-HH. The filled curves and 
dashed curves exhibit the averaged temporal profiles and their AC traces, 
respectively. The yellow highlighted region shows the tinne interval 
measured in the experiment. 



peak intensity of the main pulse reaches its maximum height at a 
relative phase of 0.1471, and the pre- and post-pulses are most 
effectively suppressed at this symmetrical point with peak heights 
of both 0.08. We find that the highest peak of the pre- and post- 
pulses in the ± a region around <(/>o) is lower than l/e^~0.14, 
whereas the area of the normal distribution in our 0o range, 
which accounts for 68.3% of all events, is depicted as grey-filled 
area in Fig. 4b. In the second-order AC experiment, the AC 
intensity of main pulse (Iac - main) and the side pulse (Iac - side) 
are given by the square of the main pulse of the CEP -averaged 
temporal profile (4^^^) and /main(^pre + ^post)> respectively. Our 
simulation shows that the event rate of the ± a region, which has 
weak pre/post pulse in Fig. 4b, contributes ~ 85% to the whole 
AC trace. From this evaluation, we can conclude that the AC 
signal at the foot of our relative phase distribution (outside the 
± a region) is negligible for composing the AC trace. In other 
words, the two -photon nonlinear process automatically picks up 
the most intense lAP part for specific combinations of the CEP 
and relative phase. 

As our method has the advantage that the HH output yield can 
be linearly scaled up by increasing the HH emission volume, we 
now discuss a scaled-up configuration in the soft-X-ray region to 
obtain even shorter lAP durations. In the soft-X-ray region 
at around 100 eV, we have already demonstrated an ~30nJ 
energy per harmonic order from OC-HH in Ne gas^^. By 
straightforwardly upgrading the TC-HHG to a main pump 
energy of 50 mj (with = 0.15) and adopting a focusing length of 
5 m with a 5 -cm-long Ne medium, we can expect to achieve HH 
pulse energies >0.1jiJ from 95 to llOeV, which is almost 
1,000-fold higher than the energies previously reported^ ^ Our 
simulations indicate that lAPs can be generated with sub-300 as 
duration. In addition, a straightforward upscaling of the lAP 
power to a few tens of gigawatts is feasible by extending the 
focusing length to 50 m with a 100-mJ main pump energy (with 
(^ = 0.15), while maintaining a pulse duration of 500 as at 
-30eV. 

Since the first demonstration of an lAP was reported in 2001 
(ref. 1), continuous effort to increase the lAP energy^^ has been 
exerted by many research groups around the world over the past 
decade, as illustrated in Fig. 7. Unfortunately, the output energy 
of lAPs has not exceeded the nanojoule level thus far. Owing to 
the low photon flux, the temporal characterization of the lAP was 



10-^ 



10- 



^ 10-^ 

CD 
CL 
< 

- 10"^ 
Q. 

° 10-^ 



10 





Quasi-IAP; 20 eV, 375 as A 




IR-TC (30 fs, w/o CEPS) 




lAP; 30 eV, 500 as 0 


28 eV, 950 as [23] 

IG (8 fs, w/o CEPS) 0. 


38 eV [37] / 

DOG (9 fs, w/ CEPS) ^y27 eV, 1 50 as [1 3] 


90 eV, 650^5 fl] 
OC (7 fs,^/oCEPS) 


80 eV, 80 as [11] 
OC (3 fs, w/ CEPS) 

35 eV 1 30 as [36] ^ ^ 40 eV, 260 as [38] 
PG(5fs,w/CEp4 • GDOG(20fs,w/CEPS) 
• 



2000 2002 2004 2006 2008 2010 2012 
Year 

Figure 7 | Progress of the energy of lAPs used in attosecond 
experiments. Red circles: previous milestone results using various 
schennes; blue points: our HHG source. IG, ionization gating; PG, 
polarization gating; DOG, double optical gating; GDOG, generalized DOG. 
We indicate the required laser pulse duration and CEP stabilization 
(CEPS) for each generation scheme. Note the dramatic leap in lAP energy, 
which is used in an attosecond experiment, by two orders of magnitude 
achieved in this work. 
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generally achieved using the attosecond streaking method 
employing frequency- resolved optical gating for complete recon- 
struction of attosecond bursts (FROG- CRAB) reconstruction. We 
have now improved the available lAP energy for nonlinear optics 
experiments by two orders of magnitude compared with the 
highest output energy reported previously^^. Our method features 
good energy scalability and high conversion efficiency owing to 
the phase -matching technique. We directly characterized the 
temporal profile of our lAP by an AC method using the nonlinear 
interaction of N2. In fact, the AC traces clearly demonstrate that 
our HH source is capable of not only inducing nonlinear 
phenomena but also having attosecond temporal resolution. This 
experiment also demonstrates the capability of our technique to 
observe nonlinear attosecond dynamics in real time with lAPs by 
time-resolved spectroscopy employing a non- sequential two- 
XUV-photon process in a molecule. Moreover, our intense lAPs 
are useful as the external injector for a seeded FEL^^ in the XUV 
region. By combining our attosecond technology and seeded FEL 
technology, we might obtain an altogether novel XUV light 
source in the near future. In addition, the quasi- 1 AP from the 
mid-plateau of TC-HH will be useful for several applications, for 
example, for the control of electron localization^"* by attosecond 
pulse trains, the manipulation of the HHG process^^'^^ and 
pump-probe experiments involving a few-cycle laser pulse. Now, 
the lAP energy has reached the microjoule level, which is 10,000 
times higher than the energy of the first demonstration. We think 
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Figure 8 | HH yield as a function of the wavelength ratio K^k-Jko. 

The relative intensity ratio ^ = /iR/Zsoonm is fixed to be 0.15. Dashed line: 
a = 0.7 scaling. Inset: scaling exponent a (from HH yield oc K~^) as a 
function of ^. 



that the presented method paves the way for the intense 
attosecond science frontier, which will certainly bring a great 
leap forward in research on ultrafast phenomena and nonlinear 
optics. 



Methods 

Generation configuration for intense attosecond pulses. From the point of view 
of energy scaling for lAPs, the following essential factors have to be considered: 
(i) scalability of the pump laser energy, (ii) high conversion efficiency of the HHG 
and (iii) scalability of the generation geometry of HHG. As proposed and 
demonstrated in ref. 14, our infrared-TC scheme enables us not only to relax the 
requirements for the pump pulse duration ( ~ 30 fs) and the CEP stabilization but 
also to markedly minimize the ionization of the harmonic medium. This is a major 
advantage for efficiently generating intense attosecond pulses, because neutral media 
allow us to use the phase-matching technique. We strictly considered and designed 
the configuration of our infrared-TC scheme for scaling up the HH output energy. 
The supplemental infrared wavelength is optimized to 1,300 nm to effectively 
suppress satellite pulses on both sides of the central pulse. To extract the HH energy 
efficiently, the harmonic medium should be extended to satisfy the absorption- 
limited condition^ ^ under the phase-matching condition. In the infrared-TC 
scheme, the difference in group delay between the 800-nm (io) and 1,300-nm (li) 
pulses even in Xe gas is evaluated to be just 9 as in the 12-cm-long medium of Xe 
(1 Torr). As this amount is negligible for the synthesized TC field, we can perfectly 
maintain the waveform in the harmonic medium during propagation. 

Next, to achieve the highest possible conversion efficiency of HHG in the 
infrared-TC scheme, we study the effect of ( = /iR/Isoonm) on the conversion 
efficiency by numerical simulations. The HH yield generally decreases with 
increasing driving wavelength, which is proportional to oc with 5 < a < 6 
(ref. 27). To investigate the scaling exponent a in the TC-HHG, we computed the 
HH yield by solving the time-dependent Schrodinger equation (TDSE), / ^^q^^ = 
[—\+ ^(0 +2:£(f)]iA(r, t), where V{r) is the atomic potential and atomic 
units are used. We assume that the laser pulse is polarized along the z direction. 
The TDSE is solved in spherical coordinates, so z = rcosO. We consider a hydrogen 
target and V{r) = — 1/r. The linearly polarized electric field is expressed by 



ETc=Eoexp [-2ln2(^)T cos(ft;of + 0ce) +^1 exp[ - 2ln2(^)'] 
cos(^ + + ^o), where K is the frequency ratio between the main field (coq) 
and the supplemental infrared field {K = <jOo/<jOi, K>1). The subscripts 0 and 1 
denote the main fundamental laser and supplemental infrared pulse, respectively. 
E is the electric field; To,i, (/)ce and 00 denote the pulse durations, CEP and relative 
phase, respectively. Note that the jitter in temporal delay 3ti and in relative 
phase 00 introduced by path-length fluctuations of the interferometer are 
connected via 3ti = KijyJcDQ. Here the pulse durations of the fundamental and 
infrared pulses are both 30 fs, and the total peak intensity is kept constant at 
/^(. = 6 X 10^3 W cm . To determine a for each ^, the central wavelength of the 
main pulse is kept at 800 nm and that of the infrared pulse is varied from 900 nm to 
2,500 nm. ^ is varied from 0 to 00 , the two limits corresponding to OC-HH 
with /lo and OC-HH with ii, respectively. The TDSE is solved using the method 
described in ref 28. The results shown in the inset of Fig. 8 indicate that the scaling 
exponent a rapidly increases when increasing the relative intensity ^ and finally 
approaches ~ 5.5. A high ^ is disadvantageous for the conversion efficiency, even if 
a good suppression of the shoulder pulses on both sides of the central pulse 
as well as an extension of the HH cutoff can be obtained. We investigated an upper 
limit of for which both a relatively high conversion efficiency and a good 
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Figure 9 | Schematic figure of the HHG setup and the second-order autocorrelator. The AC setup (bottom part) is composed of a pair of Si harmonic 
separators and controls the delay (At) between the two replicas of the reflected pulse. The two replica pulses with At delay time are focused by a concave 
mirror (SiC bulk or Sc/Si multilayer) onto a molecular beam of N2. The fragment ions from N2 molecules, which are introduced through a skimmer 
from a pulsed gas jet, are detected by a TOF ion mass spectrograph, which is constructed from three electrodes, a flight tube, and a microchannel 
plate (MCP). 
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suppression of shoulder pulses are simultaneously realized, and we decided to use 
(^~0.15. For = 0.15, we find that the HH yield obeys a scaling law of ~K~°-^ 
with the wavelength ratio K = /li//lo (see Fig. 8). Moreover, we emphasize that our 
infrared-TC method can universally be applied while maintaining the energy- 
scaling strategy of HHG, independent of the harmonic medium. 



Experimental setup. 800-nm and infrared pulses are generated starting from a 
homebuilt TW-Ti:sapphire laser system (30 fs, 150 mj, 800 nm, 10 Hz; Fig. 9). 
Infrared pulses with a 40-fs duration are produced by a white-light- seeded high- 
energy optical parametric amplifier^^. To correct for the different focal lengths of 
the 800-nm and infrared pulses, both pulses were focused using two separate 
focusing lenses (/=4.5m for 800 nm,/= 3.5 m for infrared). We discuss the 
stability of the interferometer part to produce the TC field in the Discussion 
section. The two beams are collinearly combined with a dichroic beam splitter and 
delivered into the target chamber through a thin CaF2 window. The focused beams 
pass through the interaction gas cell, which has two pinholes on the entrance and 
exit surfaces for maintaining the vacuum. The confocal parameters of the two 
beams are adjusted to be the same, matching the Gouy phase shift through the 
focus. The Xe target gas was statically filled into the interaction cell of 12 cm length. 
The HH signals were observed with a flat-field normal-incidence XUV 
spectrometer with a microchannel plate. A charge-coupled device camera detected 
the two-dimensional fluorescence from a phosphor screen placed behind the 
microchannel plate. 

To evaluate the HH output energy, we first directly measured the total pulse 
energy of OC-HH from 11th to 21st HH by using an XUV photodiode. The values 
of spectral sensitivity and quantum efficiency of the XUV photodiode were 
obtained from ref 30 and were calibrated with 266 -nm Q- switched YAG laser 
pulses in this experiment. To perfectly block the pump pulse and lower harmonic 
components (< 11th), we inserted a 200-nm-thick Al filter supported by a mesh- 
back in front of the photodiode. Although the filter transmission of Al can be 
estimated from the absorption coefficient, the Al filter transmission degrades with 
time because of oxidation^ ^ Therefore, the filter transmission for individual 
harmonics was calibrated using an XUV spectrometer. From the HH spectral 
distribution measured by the XUV spectrometer, we divided the total pulse energy 
among the six harmonic components. Thus, for OC-HH we obtain the relationship 
between the HH spectral intensity and the pulse energy. In our TC-HH 
experiments, we computed the area under the measured continuum spectrum by 
integration within a certain energy range, to evaluate the pulse energy in the cutoff 
region. When we evaluate the conversion efficiency of HH in the TC-HH case, we 
used the total pump energy (800 nm + 1,300 nm). 



Indicator of the AC measurement. We have adopted the N+ ion yield from 
fragmenting N2 molecules as the indicator of the second-order AC signal. We show 
a typical time-of-flight (TOF) mass spectrum of ions at a mass number m/z = 14 
for irradiation with the mid-plateau TC-HH field in the inset of the bottom panel 
of Fig. 2a. We can expect that the ion yields in the TOF spectrum contain the N ^ 
fragment ions from the dissociative state of (ref. 32) and the Coulomb 
explosion of N2 , in addition to doubly charged nitrogen molecular ions N2 
(ref 33) and/or N"^ ejected through the decomposition pathway of N2"^ having 
zero kinetic energy (KE). The AC function was measured by recording the side- 
peak yield at KE ~ 3 eV indicated by the inverted triangle (left side) as a function of 
the delay time At. The photon energy required for creating these KE peaks around 
3 eV is expected to be 32 eV or higher, as deduced from high-resolution 
coincidence measurements of electrons and N ^ ions^^. From the reflectivity of our 
SiC mirror we conclude that in our situation, the dominant contribution to create 
the side-peak signals at KE ~ 3 eV comes from non- sequential two-photon 
ionization of HH with energies from 12 to 21 eV. When we replaced the focusing 
mirror from SiC to Sc/Si, we observed the other KE peak of N at ~ 5 eV in the 
TOF spectrum excited by the HH cutoff. We adopted this ion yield as the indicator 
of the AC signal for the cutoff HH (~30eV), because a higher KE is expected to 
contain the signal from Coulomb explosion ofN2"^ ^N"^+N"^ (ref. 34) via a 
two-photon double ionization driven by the cutoff HH. As we could also measure 
the attosecond pulse train structure when the OC-HH cutoff was employed, we can 
conclude that the side-peak signal at 5 eV originates from a nonlinear interaction of 
the HH fields with the N2 molecules. As is shown in Figs 2 and 3, the appearance of 
the AC signal itself in the TOF mass spectra unambiguously proves that there is a 
nonlinear interaction of HH fields resulting in the production of N ions at each 
KE, as a linear interaction cannot provide a correlated signal of the pulse envelope. 

To discuss the detection sensitivity of our AC measurements, we evaluated the 
noise-to-peak intensity ratio of the AC trace. Here we defined this ratio according 
to the following procedure. First, we fit the experimental data points in the middle- 
bottom panel of Fig. 3 to a Gaussian function. We determined the amplitude of the 
AC peak (Vamp) by the fit Gaussian function trace. Next, we subtract the fit 
Gaussian function trace from the experimental data and evaluated the s.d. (Vsdev) 
as the background signal. The signal fluctuations are caused by fluctuations of the 
measured N"^ ion yield. The noise-to-peak intensity ratio is defined by T^sdev/^amp> 
resulting in 0.24. Note that we have also evaluated Vsdev by simply picking up data 
points in the ranges outside the correlation peak, which are ( — 1.15fs, — 0.85 fs) 
and (0.85 fs, 1.15 fs); again we obtained Vsdev/Kmp to be 0.21. Thus, we conclude 



that we can notice side-peak signals with amplitude higher than 0.24 times the 
main peak height in our AC measurements. Note that the signal-to-background 
ratio of all our (volume-integrated) second-order AC traces is lower than expected 
from the theory of the volume-integrated second-order AC. This lower ratio has 
the same origin as in our past attosecond pulse train measurement, namely the 
contribution of the single-photon ionization to the background^^. 

References 

1. Hentschel, M. et al. Attosecond metrology. Nature 414, 509-513 (2001). 

2. Krausz, F. & Ivanov, M. Attosecond physics. Rev. Mod. Phys. 81, 163-234 
(2009). 

3. Schafer, K. J., Yang, B., DiMauro, L. F. & Kulander, K. C. Above threshold 
ionization beyond the high harmonic cutoff. Phys. Rev. Lett. 70, 1599-1602 
(1993). 

4. Corkum, P. B. Plasma perspective on strong-field multiphoton ionization. Phys. 
Rev. Lett. 71, 1994-1997 (1993). 

5. Takahashi, E. J., Nabekawa, Y., Otsuka, T., Obara, M. & Midorikawa, K. 
Generation of highly coherent submicrojoule soft x rays by high-order 
harmonics. Phys. Rev. A 66, 021802(R) (2002). 

6. Takahashi, E. J., Nabekawa, Y. & Midorikawa, K. Generation of 10-|iJ coherent 
extreme-ultraviolet light by use of high-order harmonics. Opt. Lett. 27, 
1920-1922 (2002). 

7. Tzallas, P., Charalambidis, D., Papadogiannis, N. A., Witte, K. & Tsakiris, G. D. 
Direct observation of attosecond light bunching. Nature 426, 267-271 (2003). 

8. Nabekawa, Y. et al. Interferometric autocorrelation of an attosecond pulse train 
in the single-cycle regime. Phys. Rev. Lett. 97, 153904 (2006). 

9. Gauthier, D. et al. Single-shot femtosecond X-ray holography using extended 
references. Phys. Rev. Lett. 105, 093901 (2010). 

10. Togashi, T. et al. Intense coherent light in extreme ultraviolet region by free 
electron laser seeded with 13th harmonic of Ti: sapphire laser. Opt. Express 19, 
317-324 (2011). 

11. Goulielmakis, E. et al. Single-cycle nonlinear optics. Science 320, 1614-1617 
(2008). 

12. Sansone, G., Poletto, L. & Nisoli, M. High-energy attosecond light sources. Nat. 
Photonics 5, 655-663 (2011). 

13. Ferrari, F. et al. High-energy isolated attosecond pulses generated by above- 
saturation few-cycle fields. Nat. Photonics 4, 875-879 (2010). 

14. Takahashi, E. J., Lan, P., Miicke, O. D., Nabekawa, Y. & Midorikawa, K. 
Infrared two-color multicycle laser field synthesis for generating an intense 
attosecond pulse. Phys. Rev. Lett. 104, 233901 (2010). 

15. Constant, E. et al. Optimizing high harmonic generation in absorbing gases: 
model and experiment. Phys. Rev. Lett. 82, 1668-1671 (1999). 

16. Ayvazyan, V. et al. Generation of GW radiation pulses from a VUV free- 
electron laser operating in the femtosecond regime. Phys. Rev. Lett. 88, 104802 
(2002). 

17. Shintake, T. et al. A compact free-electron laser for generating coherent 
radiation in the extreme ultraviolet region. Nat. Photonics 2, 555-559 

(2008) . 

18. Lan, P., Takahashi, E. J. & Midorikawa, K. Optimization of infrared two-color 
multicycle field synthesis for intense-isolated-attosecond-pulse generation. 
Phys. Rev. A 82, 053413 (2010). 

19. Takahashi, E. J., Hasegawa, H., Nabekawa, Y. & Midorikawa, K. High- 
throughput, high-damage-threshold broadband beam splitter for high- 
order harmonics in the extreme-ultraviolet region. Opt. Lett. 29, 507-509 
(2004). 

20. Skantzakis, E., Tzallas, P., Kruse, J., Kalpouzos, C. & Charalambidis, D. 
Coherent continuum extreme ultraviolet radiation in the sub-lOO-nJ range 
generated by a high-power many-cycle laser field. Opt. Lett. 34, 1732-1734 

(2009) . 

21. Tzallas, P., Skantzakis, E., Nikolopoulos, L. A. A., Tsakiris, G. D. & 
Charalambidis, D. Extreme-ultraviolet pump-probe studies of one- 
femtosecond- scale electron dynamics. Nat. Phys. 7, 781-784 (2011). 

22. Takahashi, E. J., Nabekawa, Y. & Midorikawa, K. Low- divergence coherent soft 
x-ray source at 13 nm by high-order harmonics. Appl. Phys. Lett. 84, 4-6 
(2004). 

23. Sekikawa, T., Kosuge, A., Kanai, T. & Watanabe, S. Nonlinear optics in the 
extreme ultraviolet. Nature 432, 605-608 (2004). 

24. Singh, K. P. et al. Control of electron localization in deuterium molecular ions 
using an attosecond pulse train and a many-cycle infrared pulse. Phys. Rev. Lett. 
104, 023001 (2010). 

25. Schafer, K. J., Gaarde, M. B., Heinrich, A., Biegert, J. & Keller, U. Strong field 
quantum path control using attosecond pulse trains. Phys. Rev. Lett. 92, 023003 
(2004). 

26. Takahashi, E. J., Kanai, T., Ishikawa, K. L., Nabekawa, Y. & Midorikawa, K. 
Dramatic enhancement of high-order harmonic generation. Phys. Rev. Lett. 99, 
053904 (2007). 

27. Tate, J. et al. Scaling of wave-packet dynamics in an intense midinfrared field. 
Phys. Rev. Lett. 98, 013901 (2007). 



8 



NATURE COMMUNICATIONS | 4:2691 | DOI: 10.1038/ncomms3691 | www.nature.com/naturecommunications 
© 2013 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms3691 



ARTICLE 



28. Lan, P., Takahashi, E. J. & Midorikawa, K. Wavelength scaling of efficient high- 
order harmonic generation by two-color infrared laser fields. Phys. Rev. A 81, 
061802 (2010). 

29. Takahashi, E. J., Kanai, T., Nabekawa, Y. & Midorikawa, K. 10-mJ class 
femtosecond optical parametric amplifier for generating soft x-ray harmonics. 
Appl. Phys. Lett. 93, 041111 (2008). 

30. Krumrey, M., Tegeler, E., Goebel, R. & Kohler, R. Self- calibration of the same 
silicon photodiode in the visible and soft x-ray ranges. Rev. Sci. lustrum. 66, 
4736-4737 (1995). 

31. Powell, F. R., Vedder, P. M., Lindblom, J. F. & Powell, S. F. Thin film filter 
performance for extreme ultraviolet and x-ray applications. Opt. Eng. 29, 
614-624 (1990). 

32. Aoto, T. et al. Inner-valence states of and the dissociation dynamics 
studied by threshold photoelectron spectroscopy and configuration interaction 
calculation. /. Chem. Phys. 124, 234306 (2006). 

33. Wetmore, R. W. & Boyd, R. K. Theoretical investigation of the dication of 
molecular nitrogen. /. Phys. Chem. 90, 5540-5551 (1986). 

34. Sato, T. et al. Dissociative two-photon ionization of N2 in extreme ultraviolet by 
intense self-amplified spontaneous emission free electron laser light. Appl. Phys. 
Lett. 92, 154103 (2008). 

35. Nabekawa, Y. & Midorikawa, K. Interferometric autocorrelation of an 
attosecond pulse train calculated using feasible formulae. New J. Phys. 10, 
025034 (2008). 

36. Sansone, G. et al. Isolated single-cycle attosecond pulses. Science 314, 443-446 
(2006). 

37. Mashiko, H. et al. Double optical gating of high-order harmonic generation 
with carrier-envelope phase stabilized lasers. Phys. Rev. Lett. 100, 103906 
(2008). 

38. Feng, X. et al. Generation of isolated attosecond pulses with 20 to 28 
femtosecond lasers. Phys. Rev. Lett. 103, 183901 (2009). 



Acknowledgements 

This work was supported by the Ministry of Education, Culture, Sports, Science and 
Technology (MEXT) through a grant for Extreme Photonics Research; by the Japan 
Society for the Promotion of Science through a Grant-in-Aid for Scientific Research; 
by the MEXT through a Grant-in-Aid for Scientific Research for Young Scientists (B) 
No. 23760057; and through a Grant-in-Aid for Scientific Research (B) No. 25286074. 
Y.N. gratefully acknowledges the financial support by Grant-in-Aid for Scientific 
Research (A) No. 21244066. 

Author contributions 

E.J.T. was the main contributor in all aspects of this project; he conceived and designed 
the experiment. He also wrote the manuscript, which was polished by all the authors. 
E.J.T., P.L. and O.D.M. discussed the generation scheme for obtaining an intense lAP. 
P.L. performed the theoretical calculations. Y.N. contributed to the AC experiments. 
K.M. oversaw this project as a whole. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

Reprints and permission information is available online at http://www.nature.com/ 
reprintsandpermissions. 

How to cite this article: Takahashi, E. J. et al. Attosecond nonlinear optics 
using giga watt-scale isolated attosecond pulses. Nat. Commun. 4:2691 
doi: 10.1038/ncomms3691 (2013). 



(cc)©®© 



This work is licensed under a Creative Commons Attribution - 
NonCommercial-NoDerivs 3.0 Unported License. To view a copy of 



this license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 



NATURE COMMUNICATIONS | 4:2691 | DOI: 10.1038/ncomms3691 1 www.nature.com/naturecommunications 

© 2013 Macmillan Publishers Limited. All rights reserved. 



9 



